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We study the ratios = F(P -> eu e [y])/r(P -> fiv^y]) (P = ir,K) in Chiral Perturbation 

Theory to order e 2 p 4 . We complement the two- loop effective theory results with a matching calcula- 
tion of the counterterm, finding = (1 .2352 ± O.OOOf) x f(T 4 and R ( * } = (2.477±0.00f ) x f0~ 5 . 



Introduction- The ratio = T(P -> eD e [j])/T(P -> 
M^mW) = ^tK) is hclicity-suppressed in the Stan- 
dard Model (SM), due to the V — A structure of charged 
current couplings. It is therefore a sensitive probe of 
all SM extensions that induce pseudoscalar currents and 
non-universal corrections to the lepton couplings 
such as the minimal supersymmetric SM |2j. Effects 
from weak-scale new physics are expected in the range 
(AR e / / R e / p ~ 10~ 4 — 10~ 2 and there is a realis- 
tic chance to detect or constrain them because: (i) 
ongoing experimental searches plan to reach a frac- 

< 5 x 10" 4 0] and 



tional uncertainty of (AR^l )/R^9 

(AR^)/R[^ < 3 x 1CT 3 [4], which represent respec- 
tively a factor of 5 and 10 improvement over current 
errors (ii) The SM theoretical uncertainty can be 
pushed below this level, since to a first approximation the 
strong interaction dynamics cancels out in the ratio R e /^ 
and hadronic structure dependence appears only through 
electroweak corrections. Indeed, the most recent theoret- 
ical predictions read R { j9 = (1.2352 ± 0.0005) x 10" 4 @, 



R 



(1.2354±0.0002)x 10" 



0,and R e/ii 



{K) - (2.472± 
provide a gen- 



0.001) x 10~ 5 Q. The authors of Ref. 
eral parameterization of the hadronic effects and estimate 
the induced uncertainty via dimensional analysis. On the 
other hand, in Ref. the hadronic component is calcu- 
lated by modeling the low- and intermediate-momentum 
region of the loops involving virtual photons. 

With the aim to improve the existing theoretical sta- 
tus, we have analyzed R e /^ within Chiral Perturbation 
Theory (ChPT), the low-energy effective field theory 



(EFT) of QCD. The key feature of this framework is that 
it provides a controlled expansion of the amplitudes in 
terms of the masses of pseudoscalar mesons and charged 
leptons (p ~ m w ,K,e/A- x , with A x ~ farF* ~ 1.2 GeV), 
and the electromagnetic coupling (e). Electromagnetic 
corrections to (semi)-leptonic decays of K and 7r have 
been worked out to 0(e 2 p 2 ) @, @], but had never been 
pushed to 0(e 2 p 4 ), as required for R e / fl . In this letter we 
report the results of our analysis of i? e / M to Q(e 2 p 4 ), de- 
ferring the full details to a separate publication [IOJ. To 
the order we work, R e /^ features both model independent 
double chiral logarithms (previously neglected) and an a 
priori unknown low-energy coupling (LEC), which we es- 
timate by means of a matching calculation in large- Nc 
QCD. The inclusion of both effects allows us to further 
reduce the theoretical uncertainty and to put its estimate 
on more solid ground. 

Within the chiral power counting, i? e / M is written as: 
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The leading electromagnetic correction A 
sponds to the point-like approximation for pion and kaon, 
and its expression is well known 0, [TlT | . Neglecting terms 
of order (m e /m p ) 2 , the most general parameterization of 
the NLO ChPT contribution can be written in the form 



A (P) 



a m 



log -f + 4 F) + cf ] (m> P ) + 



a m 2 p _(p) rn 2 



log- 



(3) 



which highlights the dependence on lepton masses. The 

(p) 

dimensionless constants c 2 3 do not depend on the lep- 
ton mass but depend logarithmically on hadronic masses, 
while c^^m^/mp) — ► as m M — > 0. (Note that our 3 
do not coincide with 62,3 of Ref. because their C3 is 



not constrained to be m^-independent.) Finally, depend- 
ing on the treatment of real photon emission, one has to 
include in R e /^ terms arising from the structure depen- 
dent contribution to P — > eP e j that are formally of 
0(e 2 p e ), but are not helicity suppressed and behave as 
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(i) two-loop graphs with vertices from the lowest order 
effective lagrangian (0(p 2 )); (ii) one-loop graphs with 
one insertion from the NLO lagrangian [13( (0(p 4 )); (iii) 
tree-level diagrams with insertion of a local counterterm 
of 0(e 2 p 4 ). In Fig. [1] we show all the relevant one- and 
two- loop 1PI topologies contributing to R e /^. Note that 
all diagrams in which the virtual photon does not connect 
to the charged lepton line have a trivial dependence on 
the lepton mass and drop when taking the ratio of e and 
fi rates. We work in Feynman gauge and use dimensional 
regularization to deal with ultraviolet (UV) divergences. 





FIG. 1: One- and two-loop 1PI topologies contributing to 
R e /n to order e 2 p 4 . Dashed lines represent pseudoscalar 
mesons, solid lines fermions an wavy lines photons. Shaded 
squares indicate vertices from the 0(p 4 ) effective lagrangian. 
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a/ir (mp/m p ) 4 (mp/m e ) 2 . 
The calculation - In order to calculate the various coef- 
ficients cf ] within ChPT to 0(e 2 p 4 ), one has to consider 



By suitably grouping the 1PI graphs of Fig. [T] with 
external leg corrections, it is possible to show [lfj that 
the effect of the 0(e 2 p 4 ) diagrams amounts to: (i) 
a renormalization of the meson mass mp and decay 
constant Fp in the one-loop result A^ 2 ; (ii) a gen- 
uine shift to the invariant amplitude Tg = T(P + (p) — » 
£ + {pi) v t{Pv))- This correction can be expressed as the 
convolution of a known kernel with the vertex func- 
tion V = l/(V2F)fdx e iqx+iWy (Q\T(jJ? M (x) (V u - 
A v){y)U + (p)) (with V^Ap) = u7 M (7 5 )cO, once the Born 
term has been subtracted from the latter. Explicitly, in 
the case of pion decay one has (W — p — q, e i23 = +1) 
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(27r) d [q 2 -2q-p e + ie] [q 2 - m 2 + i, 
iV^q 2 , W 2 ) e^PqaPf, - A l {q 2 ,W 2 ) (q ■ pg» v - p»q v ) - (A 2 (q 2 , W 2 ) - A^q 2 , W 2 )) (q 2 g^ - q»q v ) 
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2p ■ q — q 2 q 2 
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To the order we work, the form factors Vi(q 2 ,W 2 ), 
Ai(q 2 ,W 2 ) and F^ (q 2 ) have to be evaluated to 0(p 4 ) in 
ChPT in d-dimensions. Their expressions are well known 
for d — 4 and have been generalized to any d [Io| . 
So the relevant 0(e 2 p 4 ) amplitude is obtained by calcu- 
lating a set of one-loop diagrams with effective local (Vi 
and Ax) and non-local (A 2 and Fy*) 0(p 4 ) vertices. The 
final result can be expressed in terms of one-dimensional 
integrals [10 ]. 

While c 2 4 and c 2 are parameter-free predictions of 
ChPT (they depend only on m^K, F n , and the LECs 
Lg^io determined in other processes [HI]), c^ contains an 
ultraviolet (UV) divergence, indicating the need to intro- 
duce in the effective theory a local operator of 0(e 2 p 4 ), 
with an associated LEC. The physical origin of the UV 

2 4 

divergence is clear: when calculating ST^ p in the EFT 
approach, we use the 0(p 4 ) ChPT representation of the 
form factors appearing in Eq. [5] (7^„ — > 7^, hPT ). While 
this representation is valid at scales below m p (and gener- 



ates the correct single- and double-logs upon integration 
in d d q) it leads to the incorrect UV behavior of the inte- 
grand in Eq. [51 which is instead dictated by the Operator 
Product Expansion (OPE) for the (VVP) and (VAP) 
correlators. So in order to estimate the finite local con- 
tribution (dominated by the UV region) we need a QCD 
representation of the correlators valid for momenta be- 
yond the chiral regime (7^„ — > 7^ CD ) . This program is 
feasible only within an approximation scheme to QCD. 
We have used a truncated version of large- Nc QCD, in 
which the correlators are approximated by meromorphic 
functions, representing the exchange of a finite number 
of narrow resonances, whose couplings are fixed by re- 
quiring that the vertex functions (7r|UA|0) and (7r|UU|0) 
obey the leading and next-to-leading OPE behavior at 
large q [14| . This procedure allows us to obtain a simple 
analytic form for the local coupling (see Eq. [TOf . 

( p) ~( p) 

Results - The results for c 2 3 4 and c 2 depend on the 
definition of the inclusive rate T(P — > lve["f])- The ra- 
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diative amplitude is the sum of the inner bremsstrahhmg 
component (Tib) of 0(ep) and a structure dependent 
component (Tsd) of 0(ep 3 ) [l2j]. The experimental defi- 
nition of R ^i^ is fully inclusive on the radiative mode, so 

that A^2p4 receives a contribution from the interference 
of Tib and Tsd, and one also has to include the effect 
of A^2p6 \Tsd\ 2 - The usual experimental definition of 



e/l corres P on d s to including the effect of Tib in 
(dominated by soft photons) and excluding altogether the 
effect of Tsd- consequently c£' ^= c[^\ 

Results for R^ - Defining L 9 = (An) 2 L r 9 (fi) , t P = 
log(mp//i 2 ) (ji is the chiral renormalization scale), 7 = 
Ai(0,0)/Vi (0,0), zt = (me/m^) 2 , we find: 
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log -f + f 2 + 2 «W - - 7 J log -£ + (0) 
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4(1 - Zi) + (9 - 5^) log + 2 7 1 - ^ + zt log 2^ 



3 J 2(1 - zt) 



log; 
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where «M is related to the 0(p ) pion charge radius by 

1 (AnF) 2 



r.00 



4L fl 
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The function K^ n \zt), whose expression will be given in 
Ref. , does not contain any large logarithms and gives 



a small fractional contribution to c 



(tt) 
3,4- 



As anticipated, is a parameter-free prediction of 
ChPT. Moreover, we find c 2 — 0, as expected due to a 
cancellation of real- and virtual-photon effects [lij]. Fi- 
nally, encodes calculable chiral corrections (as does 
C4,(me)) and a local counterterm c^ T (/i), for which our 
matching procedure 10] gives (za = m ai /m p ): 



= 3 CT (m) 



19 m 2 , 
~9(AnF) 2 

37 - 314 



Am 2 



7- 



114 



13(4ttF) 2 
174 - 114 
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log' 
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36^(1 - 4 
T ■ H~ Zy 

34 (-1 + 4) 3 



2 



log Z A 
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Numerically, using z,4 = a/2, we find c^ T (m p ) = —1.61, 
implying that the counterterm induces a sub-leading cor- 
rection to C3 (see Table fl}. The scale dependence of 
c§ (fi) partially cancels the scale dependence of the chi- 
ral loops (our procedure captures all the " single- log" scale 
dependence). Taking a very conservative attitude we as- 
sign to C3 an uncertainty equal to 100% of the local contri- 
bution ( I AC3 1 ~ 1.6) plus the effect of residual renormal- 



ization scale dependence, obtained by varying the scale 
(j, in the range 0.5 — > 1 GeV (| AC3 1 ~ 0.7), leading to 
Ac, ' = ±2.3. Full numerical values of C2 3 4 are re- 
ported in Table HI with uncertainties due to matching 
procedure and input parameters (Lg and 7 [l6|). 

As a check on our calculation, we have verified that if 
we neglect dj and pure two-loop effects, and if we use 
Lg = F 2 /(2m 2 ) (vector meson dominance), our results 



for 



» 

'2,3,4 



are fully consistent with previous analyses of the 



leading structure dependent corrections based on current 
algebra @, Il7j . Moreover, our numerical value of A^^ 
reported in Table [TT] is very close to the corresponding 



result in Ref. 
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(0.054 ±0.044) x 10" 
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For completeness we report here the contribution to 

(it) 

e 2 p e induced by structure dependent radiation: 
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Results for R , - In this case we have: 

e/n 



lm 2 p (r 2 )^ 



(AnF) 2 



-(1-7)-^ 
3 V 11 (AnF) 2 

is the 0(p 4 ) kaon charge radius. 
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where 

obtained from c^ 1 by replacing 31/24 — 7 



is 



-7/72 
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(P = 7r) 


(P = K) 


c (p) 
c (P) 
c (P) 

C 3 

4 P) (m M ) 




5.2±0.4 ig ±0.01 7 
-10.5±2.3 m ±0.53i, 9 
1.69±0.07i 9 


(7.84 ± 0.07 7 ) x 10" 2 
4.3±0.4 ig ±0.01 7 
-4.73 ± 2.3 m ± 0.28l 9 
0.22±0.01l 9 



TABLE I: Numerical values of the coefficients cif 1 of Eq. [3] 
(P — ir, K). The uncertainties correspond to the input values 
Ll{fi = m p ) = (6.9 ± 0.7) x 10~ 3 , 7 = 0.465 ± 0.005 [l|], and 
to the matching procedure (m), affecting only C3 . 





(P = 7r) 


(P = K) 


A% (%) 

AZ (%) 
(%) 
A M (%) 


-3.929 
0.053 ±0.011 
0.073 
0.055 


-3.786 
0.135 ±0.011 

0.055 



TABLE II: Numerical summary of various electroweak cor- 
rections to R, K '. 

13/9 7, by dropping the term proportional to log m^/m 2 , 
and by inter-changing everywhere else the label 7r with 
K (masses, t n — > £k , etc.). is obtained from by 
keeping only the second line of Eq. [H] and inter-changing 
the labels tt and ^ftT. The numerical values of and 

£2 are reported in Table HI 

Resumming leading logarithms - At the level of un- 
certainty considered, one needs to include higher or- 
der long distance corrections to the leading contribution 
A e 2 p 2 ~ —3a/ir log m p /m e ~ —3.7%. The leading log- 
arithms can be summed via the renormalization group 
and their effect amounts to multiplying R^J^ by Q 

/ \ 9/ 2 

(l-|flog^) 

1 + A LL = V 3 J = 1.00055 . (14) 
1 — — log — - 

Conclusions - In Table |TT] we summarize the various 
corrections to R^, , which lead to our final results: 

R^) = (1.2352 ± 0.0001) x 10" 4 (15) 
r(JP = (2.477 ± 0.001) x 10" 5 . (16) 

In the case of R, we have inflated the nominal un- 
certainty arising from matching by a factor of four, to 
account for higher order chiral corrections of expected 
size A e 2 p 4 x m 2 K / '(AttF) 2 . Our results have to be com- 
pared with the ones of Refs. @ and 0] reported in the 
introduction. While R , is in good agreement with both 



previous results, there is a discrepancy in R { , ' that goes 
well outside the estimated theoretical uncertainties. We 
have traced back this difference to the following problems 
in Ref. 0: (i) the leading log correction All is included 
with the wrong sign (this accounts for half of the discrep- 
ancy); (ii) the NLO virtual correction Aj , = 0.058% is 
not reliable because the hadronic form factors modeled in 
Ref. 0] do not satisfy the QCD short-distance behavior. 

In conclusion, by performing the first ever ChPT cal- 
culation to 0(e 2 p 4 ), we have improved the reliability of 
both the central value and the uncertainty of the ratios 
R^,' K \ Our final result for R, is consistent with the 

e/ n e/ n 

previous literature, while we find a discrepancy in R^f}, 
which we have traced back to inconsistencies in the analy- 
sis of Ref. |7| . Our results provide a clean basis to detect 
or constrain non-standard physics in these channels by 
comparison with upcoming measurements. 
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